ISAS/EP/94/152 



Correlation Functions Along a Massless Flow 



G. Delfino, G. Mussardo and P. Simonetti 

International School for Advanced Studies, 
and 

Istituto Nazionale di Fisica Nucleare 
34014 Trieste, Italy 



Abstract 

A non-perturbative method based on the Form Factor bootstrap approach is pro- 
posed for the analysis of correlation functions of 2-D massless integrable theories and 
applied to the massless flow between the Tricritical and the Critical Ising Models. 



1. The 2-D integrable Quantum Field Theories (QFT) associated to the massless Renor- 
mahzation Group (RG) flows between two different fixed points describe important phys- 
ical situations and present several reasons of interest. Relevant examples are provided by 
the 0(3) non-linear sigma model with the 6 = n topological term by the massless 
flows between two consecutive conformal minimal models and by the 2-D formulation 
of the Kondo problem P]. One of the important features of these QFT (in addition to 
their integrability) is the absence of a mass gap, i.e. their correlation functions have a 
power law behaviour both in the ultraviolet and in the infrared regions (with different 
critical exponents, though), separated by a non trivial crossover in between. Usual pertur- 
bative methods may fail to capture the physical content of this kind of theories and it is 
therefore useful to develop some alternative approach to their analysis. Such an approach 
does exist for the massive QFT, where the powerful method of the Form Factors (FF) 
allows us to enlighten their dynamics and provides, in particular, fast convergent series 
for the correlation functions [H, ^, |^ . In the light of these achievements, it is an important 
question to see whether or not the FF approach can be also successfully applied to the 
massless QFT. There are, however, several aspects that need a careful investigation, in 
particular the analytic structure of the FF and the convergent properties of the spectral 
representation series. The purpose of this letter is to study the simplest massless non- 
scaling invariant theory associated to the flow from the Tricritical Ising Model (TIM) to 
the Critical Ising Model. A suitable extension of this approach may be useful to inves- 
tigate a concrete physical problem, as for instance the computation of correlators of the 
Kondo problem along the flow between the two fixed points. 

2. The massless scattering theory has been developed in |T], |^, In a massless 2-D 
theory, the excitations consist in a set of left (L) and right (R) movers. These are the 
lowest energy states which propagate all along the RG flow and constitute the relevant 
degrees of freedom of the problem. Higher massive states are regarded as decoupled, their 
influence being eventually seen in the analytic structure of the physical amplitudes of the 
massless particles. The simplest way of approaching the massless integrable scattering 
theory is probably to think of it as a particular limit of a factorizable massive scattering 
theory. To avoid inessential complications, let us consider a massive integrable model 
whose spectrum consists of a single self-conjugate particle A of mass m which, in a suitable 
limit, gives rise to a massless theory with right-mover and left-mover particles. A technical 
important point is the different role played by the rapidity variable in the massive and 
in the massless theories. In the massive theory, the dispersion relation may be written in 
terms of the rapidity as 

p'' = mcosh^ , p^ = msinh6' . (1) 
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The elastic two-particle S'-matrix in the massive theories only depends on the Mandelstam 
variable s, given by s = ipiiOi) +^2(^2))^ = 2m^(l + cosh^), where 6 = 61 — 62. Use of 
the rapidity variable substantially simplifies the analytic structure of the amplitude by 
absorbing the threshold branch cut singularities in the original variable s. The physical 
sheet of the s-plane can be mapped into the strip < Im 6' < tt of the 6'-plane in such 
a way that the two-particle S'-matrix becomes a meromorphic function of the rapidity 
difference. In a massless theory, there are two possible scattering configurations. The 
first of them consists in the RL scattering, the second in the RR (LL) scattering. In the 
RL case, the Mandelstam variable s is the only independent relativistic invariant quantity. 
Since the mass gap is zero, the unitarity and crossing cuts of the S'-matrix in the s-plane 
join at the origin, so that the Riemann surface splits into two distinct parts: the "upper" 
( "lower" ) part consists in the half of the physical (unphysical) sheet with Im s > and 
in the half of the unphysical (physical) sheet with Ims < 0. Some care is required to 
introduce a parameterization in terms of rapidities since we need to send m to zero in 
eq. d^) but keeping the energy and the momentum finite. For this aim, it is sufficient to 
replace 6 hj (3 ± f3o/2 and to take the limits m ^ 0, (3q +00 in such a way that the 
mass parameter M = me^"/^ remains finite. According to the sign in front of Po, we have 

for right-movers , (2) 
for left-movers . 

In short, it is useful to regard a right-mover v4/j(/5) or a left-mover Al{i3) as formally 
defined by the limits Ar^l^P) = lim/jg^+oo A{P±Po/2). The Mandelstam variable s is now 
given by s = M^e^i"^^. Due to the aforementioned splitting of the s-plane, in the massless 
case we need two meromorphic functions of the rapidity difference in order to represent a 
function defined on the two sheets of the Riemann surface. Let Srl{i3) (S'_rl(/3)) be the 
values of the scattering amplitude on the "upper" ("lower") sheet above defined. Srl{(3) 
{.Srl^P)) is obtained from the value S{6) {S{—9)) of the massive amplitude as the limit 
AiOi) ^ Ar{(3i) and AiO^) ^ Al^, i.e. 

Snm=\im S(/3 + /3o) , S«l(/3) = lim S(-/3-/5o). 

The unitarity and crossing relations can be now expressed as Srl^/S) Sjil{P) = 1 and 
Srl{P) = Srl{P + ivr). Let us now turn our attention to the RR (LL) scattering. Since 
in these channels the variable s is identically zero, all the usual analyticity arguments 
of the S'-matrix theory cannot apply, a circumstance which clearly refiects the lacking of 
intuitive understanding of the massless scattering process in (1 + 1) dimensions. In spite 
of this fact, the amplitudes Srr and Sll can be formally defined in terms of the massless 
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limit of the massive theory. However, the rapidity shifts in the massless hmit cancel in 
these cases and therefore we simply have Srr{(3) = Sll{(3) = S{P). In view of this 
identity, the amplitudes Srr and Sll are expected to fulfill (in the rapidity variable) the 
equations valid for the massive theory. As in the massive case, it is useful to encode all the 
scattering information of the massless theories into the Faddev-Zamolodchikov algebra 

where = R,L, Ar{P) and Al{P) are the creation operators for the R and L movers 
respectively and Slr{Pi - P2) = Srl{P2 - A)- 

3. Massless form factors are defined to be the matrix elements of local operators between 
asymptotic states. To determine their analytic properties, let us firstly consider the 
matrix element of a scalar operator 0{x) between the vacuum and the two-particle state 
I A/i;(/?i)y4i(/?2) >• Let FjiL^P) and Frl^P) be the functions which, for real values of 
P = Pi — P2, take values on the upper and lower edges of the unitarity cut, respectively. 
They are related each to the other through the equations 

FrlW) = SRL{f3)FRL{P) ; (3) 

FRL{(3 + t7l) = FRL{(3-t7l) . 

Using now Flr{(3) = Frl{—[3), and the fact that Frr{(3) and Fll{P) share the same 
properties of the FF of the massive theory, the monodromy equations for the massless 
two-particle FF may be written in a compact form as 

F^,^,{(3 + 2m) = F^,^,{-(3) , 
where = i?, L. These equations are easily generalized to the n-particle FF, 

Poii...aiai^\...an iPli • • • i Pii Pi+li • • • 5 Pn) 

= SoiOi+iiPi A+l)-^oi...ai+iOi...«n (/^l) • • • 5 A) • • • 5 /^n) ; (5) 

Faia2...an{Pl + 27ri, P2, ■ ■ ■ , Pn) = Fa2...a„aA(^2-, ■ ■ ■ , Pn, Pi) , 

where . . . , =< 0\O{0)\Aaj^{Pi), . . . , Ao,„{Pn) > are meromorphic functions 

of rapidities defined in the strips < Im /3j < 27r. In the massive case, the form factors 
present simple pole singularities associated either to the bound states in the scattering am- 
phtudes or to the particle-antiparticle annihilation processes However, stable bound 
states are usually forbidden in massless theories due to the absence of thresholds. This 
leads to exclude the presence of the first kind of singularities in the RL sub-channels, which 
are the only ones to which standard S'-matrix theory can be applied. On the other hand. 
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the RR and LL scattering amplitudes formally behave as in the massive case. Hence, it 
is natural to assume that, whenever the amplitude SaaiP) has a simple pole for (3 = iu 
{u G (0, tt)) with residue g, this induces a corresponding equation for the FF 

With reference to the annihilation poles, they may only occur in the RR or LL sub- 
channels. The residue equation in this case reads 



— ireSi3/=i3+i^Faaai...a„iP', P, Pi, ■ ■ ■ , Pn) — " Yl^^^iaiPi — P) j Fai...a„iPl, ■ ■ ■ , Pn) ■ 

(6) 

The invariance of the massless theory under the spatial inversion transformation will 
provide the additional relation . . . = F'p[a„]...v[ai]{-Pn, • • • , -Pi) where 

V[R] = L, V[L] = R. 

4. A perturbation of the TIM by its sub-leading energy operator e' of conformal dimen- 
sions (3/5, 3/5) induces a massless flow that ends to the Critical Ising Model |lOl- This 



flow takes place along the self-dual line of the phase diagram analysed in |10| and a low 
energy effective Lagrangian of the corresponding QFT is given by 

= ^pd^ + ^pd^p- ^ + . . . . (7) 

This Lagrangian describes the critical Ising model (free massless Majorana fermion) per- 
turbed by the irrelevant operator {%l)dil)){il)d'4)) ~ TT, which is the lowest dimension 
nonderivative field invariant under Z2 and duality transformations. Discrete symmetries 
of the theory are ip ^ ^ —ip and ip ~^ which can be identified as 

the duality transformation and the spin reversal, respectively. The factorized massless 
scattering theory of this flow has been proposed in 0. The basic assumption is that the 
massless neutral fermions appearing in (0) are the only stable particles of the theory, i.e. 
the spectrum is the same as the infrared fixed point theory. Since Srr and Sll cannot 
vary along the flow, they are given by the commutation relations of the Ising fermionic 
fields, i.e. Siiji{P) = Sll{P) = —1, whereas for Sjil we have SiiL{P) = tanh — 
Hence, on the "upper" sheet of the s-plane we have S{s) = f^f^- Using the unitarity 
equation, this expression implies that, while the physical sheet is free of poles, the unphys- 
ical one contains two poles at s = ±iM^ which can be interpreted as resonances. Let us 
turn now to the computation of the FFQ of a local operator O. They can be parameterized 

iln the sequel, F,,i(/3i, ..,/?,; ..,/?;) =< 0\O{0MBiM ■ ■ ■ AR{f3r)AL{P[) . . . Al{P[) >. When /?i > 
P2 > ■ ■ ■ > and P'l > > ■ ■ ■ > P'l, these functions coincide with the physical matrix elements 
computed on the m-asymptotic states. 
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as 



FrAPu P2,---, Pr, P[, (3'2,---, P'l) = ^^r,/ Qr,li^l, 3^2, • ■ • , Xr] Vl, 2/2, • • • , Vl) X 

X n ^-^fT^fiiihdi3.~B',) n ^'f\z^^ . (8) 

l<i<j<r -r i=ij=i l<i<j<l yj 

where Xi = yi = e~^'' and Hr^i are normahzation constants. The auxihary functions 
III and f^L completely take into account the monodromy properties of the FF (|^). 
They are the minimal solutions of the equations faia-ziP) = Saia2{P)faia2{P + 2'7ii), with 
neither poles nor zeroes in the strip < Im/9 < 27r. Explicitly, = = sinh|, 

and 

'f3 .c^ c/tsin^ 



' {iTr-l3)t ^ 

t sinh t cosh | 



The function fuL^P) also satisfies the equation fuLiP ± iT!')fRL{P) = i±ie-t> ' where 
7 = \/2e^^^^ , G being the Catalan constant. Since 5*^^^ and Sll in this theory are 
free of poles, the FF only present kinematical poles which are explicitly inserted in (§) 
through the factors Xi + xj and i/i + yj in the denominator. With the requirement that 
the FF are power bounded in the momentum variables, Qr,/ have to be rational functions, 
separately symmetric in the {xi} and {yi} with at most poles located at Xj = or 
yi = 0. Inserting the parameterization (H) into the residue equations (P) and using 
Hr,i = —i~'^2~^^''^j'' Hr+2,i = —i~'''2~^''~^Y Hr,i+2, the recursive equations for Qr^i are given 
by 

Qr+2,/(-a;, xi, . . . , Xr] . . . , yi) = x'' {-^^f^kQrA^ii . . . , x,. ; ?/l , . . . , yi) (9) 

A ' 

Qr,i+2{.xi, . . . , x^; 1/1, . . . , yu y, -y) = y^''^^—J2 PkQrA^i, • • • , a;^ 2/1, • • • , , 

where the primed sums run over odd indices if (r + I) is even and vice versa, and pk (A;) 
are the elementary symmetric polynomials in the variables {xi} ({?/«}). 

5. Eqs. are quite general since they were obtained without reference to any particular 
operator. Here we will restrict our attention to the FF of some operators of particular 
physical relevance, namely the trace of the stress-energy tensor G(x) = T^(x), the order 
and the disorder operators $(a;) and of conformal dimensions (3/80, 3/80) in the ul- 
traviolet limit and (1/16, 1/16) in the infrared one. Since these fields are spinless, under a 

A A / r(r-l) '('~^-) ')A 

Lorentz transformation they satisfy {^ yi}) = 6 ^ ^ Qr,/({2^j}; {l/i})- 

Selection rules for the FF are obtained by assigning the transformation properties of the 
massless particles under the discrete symmetries of the theory. According to the in- 
variances of the fermionic Lagrangian (^, we assign odd Z2-parity to both right and 
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left-movers, and even (odd) parity to right (left) movers under duality transformation. 
The trace of the energy-momentum tensor is expressed in terms of the ultraviolet per- 
turbing field and the conjugate coupling constant A by the relation Q{x) = ^TTXe'{x). 
Since the subenergy e' is even under both spin reversal and duality transformation, O 
will have nonvanishing form factors Fr^i only for even r and /, starting from r = I = 2 
(the vacuum expectation value Fq^q must be identically zero since it vanishes in the in- 
frared limit). The conservation of the energy-momentum tensor implies the factorization 
Qr,i{{xi}; {Vi}) = pi\iTr^i{{xi}; {yi}). The leading infrared behaviour of ^2,2 is easily 
computed by using the Lagrangian (|^) 

F,,M,P,; P[,P',) -AnM' sinh sinh e^^+^-^i"^^ , (10) 

and, with this extra piece of information, we can completely fix its exact expression 

F2,2(/3i, (32-, (3[, (3',) = ^ sinh [] - P',) sinh . (11) 

^ ij=l,2 ^ 

The recursive equations (^ are then iteratively solved by using ([Tl|) as initial condition. 
Using H2n,2m = 7rM2i"("+i)+"*("^+i) 22("'+"^')-"-"*7-2"'^, the first right chains are deter- 
mined to be 

72„,2({a;»}; {?/»}) = (i) 

/ „ \ n-2 n-1 

T2nAM; {y^}) = (0" E P2k^ixr'-'^l , (12) 

Using now the right-left symmetry relation Tr^i{{xi}] {i/i}) = Ti^riiVi}] {xi}), one can 
immediately obtain the solution for the corresponding left chains. Let us now consider 
the magnetization operator. The invariance of the theory under spin reversal implies that 
$ has nonvanishing FF only on an odd number of particles whereas $ only on an even 
one. Taking Fi = -Fo,i = 1 as initial conditions for the recursive equations for $ and 
-^0,0 = -^1,1 = 1 for l>, we find^ 



Qr,0 








Pi"-' 




Qr,l 






Qr,2 


= p'r'^" t 'p^^t'^'^^" , 


(13) 



fc=0 



^As explained in the first reference of ||], the recursive equations for the FF of the disorder operator 
$ are given by (^) but with a plus sign in front of the product of S'-matrices. 
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Qr,3 = ^,,/2-l Yl 'Pk^-i 



where the prime denotes a sum only over even indices and r should be chosen such that 
r + / is odd for $ and even for $. 

6. The spectral representation of the two-point functions is given by 

<owo(o)>^ t ^ /;; ■ ■ ■ ■ ■ ■ "^'V ./ft. . . . . ft; . . . ft-) p 



X exp 



r,l=0 

Mr 



(14) 



where we used euclidean invariance to set x = {ir, 0). This expression clearly shows that, 
contrarily to the massive case, the convergence in the infrared limit jSi — oo,/3,- +00 
is no longer guaranteed by the exponential factor inside the integrals and completely 
relies on the behaviour of the form factors F^,/ in this limit. Let's first analyze the 2- 
point function of 0. As explicitly shown by eq. (0), -^2,2 goes exponentially to zero 
in the infrared limit so that the 4-particle contribution of Ge{r) =< 0(r)0(O) > is 
convergent. It is also easy to check that the terms with higher number of particles are 
convergent and sub-leading with respect to the 4-particle contribution. Then, plugging 
expression ( jlOD into the integral (p!^, we find for r — > 00, Gq^v) ^ n^Mir» ^ which 



perfectly matches with the expected infrared power law behaviour. On the other hand, 
the behaviour of this function in the ultraviolet limit r — is fixed by the conformal OPE, 
i.e. Ge(r) ~ (|7ra)^ (Mr) "^^/s^ where a = 0.148695516 §. A logarithmic plot of 
G@{r) obtained by including the first two contributions in the spectral representation 
(|np is given in fig. 1. The slope of the curve interpolates between the two values —8 and 
— 12/5 relative to the infrared and ultraviolet fixed points; the figure also shows a very 
fast ultraviolet convergent pattern, analogously to the massive cases. The fast convergent 
behaviour of the series is also confirmed by the c-theorem sum rule Ac = ^ J drr^ < 
0(r)0(O) > Inl. In fact, we obtain Ac^^^ = 0.19600 ± .00007 only including the 4-particle 



contribution and Ac''^^ = 0.1995 ± .0004 adding the 6-particle one, where the expected 
value is Ac = 7/10 — 1/2 = 0.2. The fast convergent pattern of the spectral series 
also occurs for the 2-point function of the energy operator e, with conformal dimensions 
(1/10, 1/10) in the ultraviolet limit and (1/2, 1/2) in the infrared one, as shown in fig. 2. 
All its FF can be computed by the recursive equations (^ with the initial condition 

A more subtle situation occurs for the correlation functions of the magnetization op- 
erators. The leading infrared behaviour of the FF of $ and $ can be determined by 
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specializing the FF bootstrap equations to the Ising critical point where Sm = 
fRLiP - 13') = e(^-^')/2. With the initial conditions Fi,o = Fo,i = Fo,o = 
limiting expressions of the FF are given by 



K,l (A) ■ ■ ■ y Pr, Pi, ■ ■ ■ 



tanh 

i<j 



tanh ■ 



-1 and 
1, the 



(15) 



where r + / is odd for $ and even for $. Expression (|T3p clearly shows that the integrals 



in the spectral representation of G<s,{r) =< $(r)$(0) > and 



--< $(r)$(0) > are 



infrared divergent, i.e. these correlators require the resummation of the whole (suitably 
regularized) spectral series. Although, in general, it is an interesting open problem to 
develop mathematical tools that allows us to deal with such series, here we want to show 
that an exact resummation can be performed in the low energy limit so that the exact 
infrared conformal dimensions A$ = A|, = 1/16 can be extracted. In fact, using the FF 
(llSD and exploiting their complete factorization in a left and a right part we can write 



G^«(t)^G4«(t)+Gf(t) 



^IR 



IR/ 
I- 




^1=1 



r=l 



where t = Mr/2 and A is an infrared cutoff. Shifting the rapidities, it can be expressed 
as 



G'^^(t) 




-te- 



(16) 



As in 0, this can be considered as the square of the partition function of a classical gas of 
particles living on a semi-infinite line and subject to the pairwise interaction V{l3i — I3j) = 



— In (^tanh 



The activity of such a gas is given by z{j3) 



2w 



This function 



has a plateau zq = l/27i inside a box of length L ~ In Hence, removing the cutoff A is 
equivalent to take the thermodynamic limit L ^ oo of a gas with a constant activity Zq. 
Standard relation between the partition function and the bulk-free energy per unit length 



f{zo) allows us to write G (t) 



~ e 



2/ 



2/^ 



The bulk-free energy f{zo) can 
exactly be computed for generic values of f^], /(^o) = ^ arcsin(27r2o) — 2~2 arcsin^(27rzo), 
so that /(l/27r) = 1/8, which coincides with the expected infrared anomalous dimension 
of the fields $(x) and $(a;). 

As explicitly shown in this paper, the FF of massless theories can be exactly computed 
and they constitute useful non-perturbative information of the RG flow. The above anal- 
ysis also shows that there are sectors in the theory where the spectral series presents 
a fast rate of convergence in the entire domain of definition while for others the power 
law behaviours are reproduced by summing the whole series of FF. We expect that this 
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pattern is also present in other massless QFT and therefore, to take full advantage of the 
non-perturbative knowledge of FF, it would be interesting to develop new resummation 
methods of the spectral series. 

The authors are grateful to F.A. Smirnov and A.B. Zamolodchikov for useful discus- 
sions. 
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Figure Captions 



Figure 1 . Logarithmic plot of G@{r) including the 4 and 6 particle contributions. 
Dashed line: leading OPE of the ultraviolet fixed point. 

Figure 2 . Logarithmic plot of the 2-point function of the energy operator e including 
the 2 and 4 particle contributions. 
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